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Abstract 

In this work we compare two different random dilutions on a mean 
field ferromagnet: the first model is built on a Bernoulli-diluted graph 
while the second lives on a Poisson-diluted one. While it is known 
that the two models have, in the thermodynamic limit, the same free 
energy, we investigate on the structural constraints that the two models 
must fulfill. We rigorously derive for each model the set of identities 
for the multi-overlaps distribution, using different methods for the two 
dilutions: constraints in the former model are obtained by studying the 
consequences of the self- averaging of the internal energy density, while 
in the latter are obtained by a stochastic-stability technique. Finally 
we prove that the identities emerging in the two models are the same, 
showing robustness of the ferromagnetic properties of diluted graphs 
with respect to the details of dilution. 

Keywords: diluted graphs, spin glasses, polynomial identities. 

1 Introduction 

In the past decade an increasing interest has been shown in statistical mechan- 
ics built on diluted graphs (see i.e. [3] |I0] [12] [27] [28] [30] ) . For diluted spin 
glasses [21] [7] this interest is, at least, double motivated: despite their mean 
field nature, they share with finite-dimensional models the fact that each spin 
interact with a finite number of other spins. Moreover they are mathematically 
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equivalent to some random optimization problems (i.e. K-SAT or X-OR-SAT 
depending on the size of the instantaneous interaction |25j|26] ). 

Although the cases of simpler ferromagnetic models P]f5]fT^ are not in- 
teresting from the point of view of the hard satisfiability interpretation, they 
are still interesting for their finite connectivity nature and for testing different 
variations on the topology of the graph they are based on. 

With this aim we consider two different ways of diluting the graph [13] : the 
first ferromagnet has its links distributed according to a Bernoulli probability 
distribution, the second according to a Poisson. 

For these models we compared the properties of a family of linear con- 
straints for the order parameters (known as Aizenman- Contucci polynomials 
[2] [6] in the case of spin glasses) . 

These relations were investigated earlier [2] [21] in the spin glass framework, 
where they were obtained as a result of the stability of the quenched measure 
with respect to random perturbation, or equivalently through the bound on 
the fluctuation of the internal energy. Here we propose them as a test for 
robustness under dilution. 

The methods to approach the identities for the two models are structurally 
different. For the Poisson case, in fact, the additive law of Poissonian ran- 
dom variables makes possible the direct exploitation of the stochastic stability 
property. The same strategy is not available for the Bernoulli random vari- 
ables but, for those, we derive the set of identities from the general bound 
on the quenched fluctuations (even though, for the sake of completeness, in 
the Appendix we derive the constraints within this general framework for the 
Poisson case too). The methods we use are generalizations of those appearing 
in [83[llj[16JLl7ji23j[29j[15j[20j[6j[r]. 

Our main result is a rigorous proof of the identities and, especially, the fact 
that they coincide for the two dilutions. 

2 The mean field diluted ferromagnet 

We introduce a large number of sites, labeled by the index i = 1, ...,N, and 
associate to each of them an Ising variable (jj = ±1. 

We introduce furthermore two families of discrete independent random vari- 
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ables uniformly distributed on 1,2, ...,N. 

The Hamiltonian Hn{cf) of the diluted ferromagnet is expressed trough 

X 

HN{a) = - ^ai^aj^ (1) 

where x does depend on the dilution probability distribution. 

For the Bernoullian dilution case the variable x is called k, and defined by 

(2) 

M — N(N — l)/2 being the maximum amount of couples a^aj existing in the 
model and a/N the probability that two spins interact. 

a > plays the role of the connectivity. 

The mean and the variance of k are obtained as 

Eslk] = (3) 
Eslk'l-Ellk] = ^(1-^). (4) 

We will make use in the sequel of following properties of a Bernoulli distribu- 
tion: 

Mkgik)] = ^EB[gik + l)], (5) 
EB[k'g{k)] = M^^_iylEB[g{k + 2)]-^EB[g{k + l)], (6) 
^Mg{k)] - ^Es[g{k+l)-g{k)]. (7) 

For the Poissonian dilution case, we denote x as ^^jv, which is a random 
variable of mean aN, for some a > (again defining the connectivity of the 
model), such that 

PiCaN = k) = 7r{k, aN) = exp ( - aN^ k^O,l, 2, ... (8) 
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Furthermore, we stress that the Poisson distribution obeys the following prop- 
erties 

k7r(k,X) = A7r(A;-l,A), (9) 
^7r(fc,A) = -7r(fc,A) + 7r(A;-l,A)(l-M- (10) 

As for the Bernoulli case, the average with respect to the Poisson measure will 
be denoted by the index 

fe=0 

The expectation with respect to all the quenched variables will be denoted by 
E and represents the product of the expectation over the dilution distribution 
and the expectation over the uniformly distributed variables 



E - Eb,p 



The thermodynamic objects we deal with are the partition function 

Z^,K/3) = 5^e-/^^-("), (12) 

the quenched intensive free energy 

^iv(«,/3) = ^ElnZjv(Q;,/3), (13) 

the Boltzmann state 

^(^(<^)) - E ^('^)e"'''~^"^ (14) 

the replicated Boltzmann state 

Q{g{a)) = l[u;^^\g{a^^^)) (15) 
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and the global average {g{cr)) defined as 

{g{a)) = E[n{g{cr))]. (16) 



The functional order parameter of the theory is the infinite family of multi- 
overlaps, defined as 

N 



1 V- ^ 



i=l 



where particular emphasis is due to the magnetization m = gi = (1/A^) J2iLi '^i 
and to the two replica overlap qi2 = (1/A^) Yli=i 



3 BernouUian diluted case 

Identities in the BernouUian model will be obtained as a consequence of the 
internal energy self-average [H] [20] ; before focusing on this procedure, let us 
recall that 

Definition 1 A quantity A(cr) is called self-averaging if 

hm {(A{a) - {Aia)))") = lim f{A\a)) - (A(a))^) = 0, (17) 
by which we can remind the 

Proposition 1 Given two regular functions A{a) and B{a), if at least one of 
them is self -averaging, then the following relation holds 

lim {A{a)B{a)) = lim {A{a)) lim (5(a)) (18) 

N^oa N—*oo N—i-oo 

Proof 

The proof is straightforward. Let us suppose the self-averaging quantity is 
B{a) and use A{(7) as a trial function. Then we have 

< \{A{a)B{a))-{A{a)){B{a))\ (19) 
= \{ Aia)Bia ) - A(a)(g(cr)) + {A ia))B{a)) - (A(a))(i?(a))| 
< ^/{^)^{iB{a)-{B{aW), 
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where, in the last passage, we used Cauchy-Schwartz relation. 
In the thermodynamic the proof becomes completed. □ 

The scheme to follow is then clear: using the above proposition as the un- 
derlying backbone in the derivation of the constraints in this section, we must, 
at first, show that the internal energy density of the model self-averages. Then 
we use as trial functions suitably chosen quantities of the order parameters. 
The identities follow by evaluating explicitly both the terms of eq. (fT8|) : this 
operation produces several contributions, all involving the order parameters, 
among which cancelations occur and the remaining part gives the identities. 



3.1 Self-averaging of the internal energy density 

Once defined hi = H{a^''^)/N as the density of the Hamiltonian evaluated on 
the generic l^^ replica and 

e = tanh(/3) (20) 
a' = Ma/N^ a/2 (21) 

for simplicity, we start with the following 

Theorem 1 In the thermodynamic limit, and in (3 -average, the internal en- 
ergy density self-averages 

lim r\{n{h'^)-n{h)Ad(3 = ^. (22) 

This Theorem has already been proved in [18] by using essentially the existence 
of the thermodynamic limit for the free energy density. To make the paper 
self contained we provide an alternative proof in the appendix. 

We can now introduce the following lemma. 
Lemma 1 Let us consider for simplicity the quantity 

s 

AG = J2 [E{n{hiG) - Q{hi)Q{G)) . (23) 

1=1 
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For every regular, smooth function G, in /3 -average, we have 



rP2 

lim / \/\G\d(3 = Q. 



(24) 



Proof 



/ \AG\dp < / XI I H^i^iG) - n{hi)n{G)] I di3 

hi hi 1=1 

r02 * 

< / V ^/E[{n{hiG) - n{hi)n{G)y]d(3 (26) 
hi 1=1 

rf32 

< s / VE[^](/^2) -^]2(/^)]^^ (27) 

Jpi 

< 8^132 -Pi J p E[fi(/i2) - n'^{h)]dl3 ^ (28) 



(25) 



where (1251) comes from triangular inequality; fl26|) is obtained via the Jensen 
inequality applied to the measure E[-]. In the same way fl27|) comes from 
Schwarz inequality applied on the measure fi(-) (being G bounded), while fl25]) 

is obtained via Jensen inequality applied on the measure {P2 — Pi)~^ 
□ 

Now we can state the main theorem for the linear constraints. 

We are going to introduce a specific trial functions that we call /^(a,/?). 



Theorem 2 Let us consider the following series of multi- overlap functions G 
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acting, in complete generality, on s replicas 



1=1 

l,s 



a<l 



l<a<b l<a 



2 ^ - -Js+l,s+2 

s{s + l) 



3! 



,8+1,5+2) + 
(29) 



in the thermodynamic limit the following generator of linear constraints holds: 



P2 

lim / d(3\fG{a,(3)\ = 0. 



(30) 



Proof 

Let us consider explicitly the quantities encoded in (!23l) . For the sake of 
clearness all the calculations are reported in appendix, here we present just 
the results. 

s 

E[n{hiG)] = -a'[{Gm^) + e[Y.{Gql)-s{Gql^,)) + 

a=l 

+(^'{i:{Gql,) - sJ^iGql^,^,) + '-^^{Gql^,,^,)) 

a<b a 

+0{9')], (31) 
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-a'[{Gmf)+e[j2{Gql) - (s + l){Gql^,)) + 

a=l 

l,s l,s 

+e'[j2(Gml) -{s + l){Gm^) + Y^{Gql,) + 

a a<b 

-{s + 1) j:{Gql,^,) + 



+0{9'')\. 



(32) 



Subtracting the last equation from the former, immediately we conclude that 



AG = -fG{a,P), 



(33) 



from which theorem thesis follows. □ 



3.2 Linear constraints for multi-overlaps 

We outline here the first order identities as it is customary to do in the spin- 
glasses counterpart [6] or in neural network [9]. 



Proposition 2 The first class of multi- overlap constraints is obtained by choos- 
ing G = . 

In fact, if we set G = qf = m^, the function fG{a,l3) becomes 



fm<a,P) = a'[((mt)-(mX))(l-e2) + 
-29^{mlql,)-{mlql,)) + 
+3e'({mlql,,)-{mlql,,))+0{e') 
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from which, changing the Jacobian dO = {1 — 6'^)d(3, we get 



hm 



\fm^{a,(3)\d(3 



Pi 



d- 



1-^2) 

-2e[{mlql,)-{mlql,)) + 



{m\) - {mimi) + (34) 



0, 



where, the (not interesting) breakdown at = 1, of the expression above, 
reflects the lack of convergence of the harmonic series we used in eq. (162!) . 



Proposition 3 The second class of multi- overlap constraints is obtained by 
choosing G = g^^. 

In fact, if we set G = qf2, the function fG{a,P) becomes 

+2e(^{qt2)-Mql2ql3) + Hquql)) + 

Again 



lim 



\fAc^,P)\dp 



d- 



K?i2)-K?r2) +(35) 



1-^2) 

-3^'((g?2g?23) - 3(g?2g234) + '2^q\2<i,^ 

+0(^3) il = 0, 



from which the constraints are obtained as the r.h.s. of fl34|35l) set to zero. □ 
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4 Poissonian diluted case 

To tackle the Poisson diluted ferromagnet, we are going to use the cavity field 
approach. 

The idea beyond this technique is that, calling F{j3) the extensive free energy 
and f{/3) the intensive one, a bridge among the two, in the large N limit, is 
offered simply by the relation 



As our system has a topologically quenched disordered the + 1 spin can be 
seen as an "external random cavity field" for the former system of particles. 
The identities derived by tuning that field are called of "stochastic stability". 
The simplest way to find them is to consider monomials which are left invariant 
by the random field: the derivative with respect to it, being zero, will produce 
the desired polynomial. 

4.1 Cavity field decompositions for the pressure density 

To start applying the sketched plan let us decompose (in distribution) a Pois- 
sonian random Hamiltonian of A^ + 1 spins in two Hamiltonians [1] : The former 
of the "inner" A^ interacting spins, the latter as the remaining spin interacting 
with the internal A^ spins of the cavity. 

Up to negligible corrections that go to zero in the thermodynamic limit, 
we can write in distribution 



F^+i(/?)-F^(/?)) =/(/?) + 0(Ar-i) 



(36) 





(37) 



or simply for compactness 



(38) 



where 



N 
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It is useful now to introduce an interpolating parameter t G [0, 1] on the term 
encoding the linear connectivity shift so to menage the derivative with respect 
to the random field by differentiating with respect to this parameter. 

Definition 2 We define the t-dependent Boltzmann state uJt as 

uM^)) = ^ ] Y.aW)^^'''''''^'^'''- (40) 

We stress the simplicity by which the t parameter switches among the system 
of + 1 spins and the one built just by the former iV, in the large N limit: 
In fact, being the two body Hamiltonian left invariant by the gauge symmetry 
ai eat for all i G (1, N) with e = ±1, by choosing e = a^^i we have 

^7v,t=i(a,/3) = ZN+i{a,p), (41) 
ZN,t=oia,/3) = ZMia,/3). (42) 

Note that Z'7v,i(a, (3) is defined accordingly to fl40p and coherently, dealing with 
the perturbed Boltzmann measure, we introduce an index t also to the global 
averages (.) {.)t. 

4.2 Stochastic stability via cavity fields 

We are now ready to attack the problem. 

We divide the ensemble of overlap monomials in two large categories: stochas- 
tically stable monomials and (as a side results) not stochastically stable ones. 
Then we find explicitly the family of the stochastically stable monomials, and 
by putting their t-derivative equal to zero we obtain the identities. To follow 
the plan let us introduce the 

Definition 3 We define as stochastically stable monomials those multi-overlap 
monomials where each replica appears an even number of times. 



We are ready to introduce the main theorem, which offers, as a straightforward 
consequence, a useful corollary, stated immediately after. 
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Theorem 3 At t = 1 the Boltzmannfaktor of the perturbed measure is com- 
parable with the canonical Boltzmannfaktor, and, in the thermodynamic limit, 
we get 

lim EQN,t=i{crh<yi2---<yin) = lim EQN+i{,cri^cri,,...ai„a'^^i). (43) 

Corollary 1 In the thermodynamic limit, the averages {■)t of the stochastically 
stable monomials become t-independent in (3 -average. 



Let us focus on the proof of Theorem [31 Corollary [T] will be produced as a 
straightforward application of Theorem [3] on stochastically stable monomials. 
Let us start the proof. Let us assume for a generic multi-overlap monomial 
the following representation 



where a labels replicas, the inner product accounts for the spins depicted by 
the index / which belong to the Boltzmann state a of the product state Q for 
the multi-overlap qa^a' and runs over the integers from 1 to the amount of times 
the replica a appears into the expression. 

The external product multiplies all the terms coming from the internal one. 
The factor / fixes replica-bond constraints. 

For example the monomial Q = q'i2Q'23 has s = 3, = n'^ = 1, ra^ = 2 and 
/ = A^~^5ji j25j2 j3, there the (5-functions give the correlations 1,2 gi_2 and 



By applying the Boltzmann and quenched-disordered expectations we get 



Proof 




2,3 



92,3- 



s 



i1 a=l 1=1 



Let us suppose now that Q is not stochastically stable (for otherwise the proof 
will be simply ended) and let us decompose it by factorizing the Boltzmann 
state uj and splitting the terms involving replicas appearing an even number 



14 



Adriano Barra, Federico Camboni, Pierluigi Contucci 



of times from the ones involving replicas appearing an odd number of times. 
Then, evaluate the whole receipt at t = 1. 

u n" s rJ' 

(g),=E5^/(K},{^n)n^«(n<?) n ^''(n^)' 

jo^ji a=l 1=1 b=u+l 1=1 

where u stands for the amount of replicas which appear an odd number of 
times inside Q. 

In this way we split the measure fl in two ensembles Ua and Ub- Rephcas 
belonging to ujb are an even number while the ones in Ua an odd number. 
At this point, as the Hamiltonian has two body interaction and consequently is 
left unchanged by the symmetry cr" — >• cr"(j^_,_;^, G (1, A^) (as crj^^i = 1), we 
apply such a symmetry globally to the whole set of N spins. The even measure 
is left unchanged by this symmetry while the odd one takes a multiplying term 

O'N+l 

u 71^ s 

ia^ib a=l 1 = 1 b=U+l 1 = 1 

The last trick is that, by noticing the arbitrariness of the N + 1 label in aN+i, 
we can change it to a generic label k for each k ^ {if} and multiply by 
1 = N~'^ Ylk=i- -^^ finite N the thesis is recovered forgetting terms 0{1/N) 
and becomes exact in the thermodynamic limit. □ 

It is straightforward to check that the effect of Theorem [3] is not felt by stochas- 
tically stable multi-overlap monomials (Corollary [1]) thanks to the dichotomy 
of the Ising spins (cr^_|_i = IVra e N). □ 

The last point missing to obtain the identities is finding a streaming equation 
to work out the derivatives with respect to the random field of the stochasti- 
cally stable monomials. To this task we introduce the following 

Proposition 4 Given as a generic function of the spins of s replicas, the 
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a=l 



a<b a=l 

+ \Fs(^io ^io /«,"J + 

l,s l,s 

a<b<c a<b 



+ 



2! 

a=l 



+ 31 {FsC^io <^io <^io )t,a\ + 0{9 ) 



Proof 

The proof follows by direct calculations 



dt at ' ^ e^i=iii3j:^z^ 'rt'Tt+pj:^if 



Y F e^»=i ^'^^o +P ^S?" < '^"^ +^ ^-'5* < ) 
2aE[ ^^"^ ' P P ] - MFs)t 

Z^{(t} 

^^.^. a(i^sn:,,(cosh/3 + < sinh/j)) ^ 
0,(nu(cosh/? + <7«sinh/3)) 



16 



Adriano Barra, Federico Camboni, Pierluigi Contucci 



then, by noting that 

1 



a=l 



a<fe<c 



;i + cui^)^ 



1 — suJtO 



S{S + 1) 2 .2 K^ + l)(^ + 2) ,..3.3 , 

-uJtU uJtU + ... 



2! 



3! 



we get 

^i^ = 2aE[a(F,(l + E<^ + E<<^'+ E <«^' + -))x 

a=l a<h a<b<c 

from which the thesis follows. □ 



4.3 Linear constraints for multi-overlaps 

We saw that the stochastically stable multi-overlap monomials become asymp- 
totically independent by the t parameter upon increasing the size of the system. 
Calling for simplicity Gn^q) a stochastically stable multi-overlap monomial, 
identities follow as a consequence of Corollary [T] and are encoded in the fol- 
lowing relation 

lim dt{GM{q))t = 0. 

As we did when we investigated the Bernoullian model, we analyze the stabil- 
ity of (m^) and (^12), up to the third order in 9, so to compare the results at 
the end. 



dt{ml)t 



+ 26^3 / (^mlqu3)t - {nilq2u)t ] + 0{9 
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+ 2a^^ Kg^23)-K?234) 



(46) 



- 4«^3j^(g?29i23)* - 3(g?2fe4)t + 4(g?2g345)*) + 



4«^' (gfsgm) - 3(gf29i34) + 2(gf2%45) + 0{e 



5 Discussion and outlook 



(47) 



Let us start this section by comparing the results we get from the two models. 
We have to compare respectively eq.s (!34|) versus (H6l) and (135!) versus (H7I) . 
We see that the details of the dilution do not affect the constraints: the series 
show the same set of identities. 

In fact, despite we are not generally allowed to set to zero each term in the 
expressions fl46|47|34ll35l) (as we do to obtain alone the following identities 



(14811531) )■ at least close to the critical line, where different multi-overlaps have 
different scaling laws [1], i.e. oc {aO — 1)", such a spreading is possible and 
we can forget each single (a, /3)-coefiicient as it does not affect the identities 
(it is never involved into the averages (.)). 
We get 



= {m\) - {Tnjml), 

= ("^1^12) - ("^??23)' 

= {mlql^^) - {mlql^^}, 
when investigating the magnetization as a trial function and 
= {mlql^) - {mlql^), 
= {qt2)-MquQl3)+^{QuQl), 

= {QuQus) - 3(g?2?D + 2(g?29345). 



(48) 
(49) 
(50) 

(51) 
(52) 
(53) 
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when investigating the two rephca overlap. 

Even if a minor point, we stress that the differences among the global {a,(5)- 
coefficients (not shown here) clearly are related to the differences in the two 
method involved (in the former the constraints appear under the integral over 
the temperature, while in the latter this /3-average is already worked out), 
furthermore the limiting connectivity in the Bernoulli dilution is a/2, while is 
2a in the Poisson model; so there is an overall factor 4 of difference among 
the results (for the sake of clearness we worked out in the appendix also the 
constraints in the Poisson diluted case via the first method, to check explicitly 
the coherence among the two model coefficients). 

Then, by looking explicitly at the constraints, several physical features can be 
recognized, in fact every term is well known: 

the first class (Eq.s H8lH9|50|) is the standard first momentum self-averaging on 
replica symmetric systems; In fact, by assuming replica equivalence, eq. fHHj) 
turns out to be the standard internal energy self-averaging of the Curie- Weiss 
model. Eq. fH9l) and fl50|) contribute as higher order internal energy self-average 
by taking into account the dilution (in fact, they go to zero whenever a oo 
because ^-powers higher than one go to zero and only the Curie-Weiss self- 
averaging for the internal energy survives as it should). 

With a glance at the identities coming from the second constraint series (Eq.s 
l5T]|52ll53p we recognize immediately the replica symmetry ansatz for the mag- 
netization in the first identity, followed by the first and the second Aizenman- 
Contucci relation for systems with quenched disorder [2] [1] . 
Interestingly these series are in agreement even with other models, apparently 
quite far away, as spin-glasses with Gaussian couphng AAfl, 1] instead of ATfO, 1] 
|15] . A very interesting conjecture may be that these constraints hold for sys- 
tems whose interaction has on average positive strength and are affected by 
quenched disorder, independently if the disorder affects the strength of the 
interaction or the topology of the interaction. We plan to report soon on this 
topics. 
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A Appendix 



A.l Alternative proof of energy self-averaging 

Starting from the thermodynamic relation 



N d(3 



we evaluate explicitly the term E\p.(}i)\ as 



E[^](/^)] 



-1e 

N 



-a'E 

-a'E 

-a'E 
-a'E 



ZN{a,P) 



{a} ^io^jo^' 



Y.{ct} o-ioO-jo(cosh/3 + Oi^^aj^, sinh/3)e" 
X;|^|(cosh/3 + o-fcoO-^o sinh/5)e-^^ 



^^(1 + ^.0^/0^)6-'^^ 



(54) 

(55) 
(56) 
(57) 

(58) 

(59) 
(60) 



.1 + uj{akoaig)9. 

where in (j56ll we fixed the index u, in (j57ll we used the property ([5]) of the 
Bernoulli distribution and we introduced two further families of random vari- 
ables {ki,},{l^}, and in (l58il we used e^'^'o'^^o = cosh/9 + (TioO"jQ sinh/5. 
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Let us now expand the denominator of (l60il keeping in mind the relation 
1 



{i + Otoy - ---- ■ 2 

such that, by posing p = 1, we obtain 



3! 



E[n{h)] = -a'E\e + J2i-ir9-{l - 9'){qlJ 



n=l 



(61) 



By applying the modulus function to the equation above we can proceed further 
with the following bound 



\E[Q{h)]\<a'E[\e\ + J2\ni-0'M...n)\ 

n=l 

Both 1^1 and belong to [0, 1] so we get 

oo 



(62) 



n=l 



(63) 



whose harmonic series converges to 1/(1 — 6), \6\ < 1; 

The fact that the convergence is not guaranteed at zero temperature with this 
technique is not a problem because, first the identities we are looking for hold 
in /3-average, secondly the zero temperature has been intensively investigated 
elsewhere |22j . 

For each finite /5, then, we can write 



\E[n{h)]\ < a' 



= a 

= a' 
< 3a' 



1 + 
1 + 



i-e 
:i-e){i + e) 



i + (i 



1-^ 



< 



(64) 

(65) 

(66) 
(67) 
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and consequently 

r02 rf^- 



[ ' E[Q{h^) - Q\h)]dp < [ ' \ E[n{h^) - n\h)] \dp = 

J Pi Jpl 

1 r^'' d 

' E[Q{h)] \df3 < 



N- 



N Jp^ ' dp 

a' , , 

< 3- (68) 



rP2 

lim / EpQv')-n'^{h)]dl3 = Q (69) 



and the proof is closed. □ 



A. 2 Details in Bernoulli dilution calculations 

Here some technical calculations concerning the self-averaging technique on 
the Bernouilli diluted graph are reported. 



E[^l{hiG)] 



1 yT< XT'" a' G e^PT.U^TX=.-t-V 

'n 



<{ .} 



jHj 

N 



and remembering the properties of the BernouUian distribution ^ we can 
write 



aM 
'1^ 



E 



E 



E{.} nLi [cosh /3 + <^^o sinh (3] e- P^^ ^ 

E{.} n:=i[cosh/5 + a-^a^ sinh/3] e'/^^^ 
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■E 



aM 



E 



1+UJ 



let us expand both the numerator and the denominator up to the second order 
in 



aM 



E 



X 



aM 



E 



X 



X 



a; 



Ar2 

l,s 



E 



aM r 



^0 Jo> 

s 



\Gm^l) + 0ij^{Gql,)-s{Gql^,)) + 



a=l 



a<6 
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While the other term 'E[Q{hi)Q{G)] can be worked out as follows: 

E[nihi)niG)] = 



N 
1 



E 



E 



7V2 



E 



aM 



E 



E 



s+1 



a a<6 



aM 



l.s 



E 



X 



a 

l,s 

a a<6 

X (|l - (s + l)uj{ai^aj^)9 + 



+ 



(s + l)(s + 2) 2. ^ ..2 



Ar2 



■E 



l,s 



(s + l)f](G'<4j(^Ka,j) + 
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a a<b 

-{s + l)n{G)u;{ai,crj,) + 

l,s 

+ 1) ^ 0(G<<7« <<7}ju;(<7,„<7,J 

a 



+ 



a=l 

( Y.^Gml) - (. + l)(Gm,2) + Y.^Gql,^) + 



a<6 



l.s 



s+l,s+2 



')] 



A. 3 Poisson identities via the self-averaging technique 

For the sake of completeness we report also the constraints in the Poisson 
diluted model obtained by using the first method: 



1=1 



+ 



a<l 

l.s 



l.s l,s 

+ 3e\j2(G<lla,,)-sY.{Gql,s+,) + 

l<a<b l<a 



s(s + 1) 



l.s 



^^^if^(G,.%,,.«,.„))+o(e')], 



(72) 
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by which, choosing as a trial function we have 



-2^(K?12)-K?23)) + 
+Se'[{mlql,,) - (mlql,,)) +0(0')], 



from which, changing the Jacobian d9 = {1 — 9'^)dp, we get 

f/32 1-62 

„ e 



hm 

AT— >oo 



ml) - {mlml)^ 

(1 _ 02^ ((^??12) - (^??23)) 



+ 
+ 



02 , 



(73) 



If we set G — (^2^ the trial function /^'(q;, /3) becomes 

+ 2^((g^2) - 4(gi22?23) + 3(g?2?34)) + 

Again 

(?^2) - 4(gL?23) + 3(gL?34)) 



lim 

AT— ►oo 



+ 



2q; / 
9 



+ 



^2 



(74) 



(1 - e^) 
+ o(^')l]=o. 



(5l25l23) - 3(^12^234) + 2(5129345) 
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